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ABSTRACT 


Electron  trajectories  are  obtained  for  the  Hahn-Ramo  space- charge 
waves.  The  approach  is  to  use  Lagrangian  form  of  hydrodynamic  equations 
of  motion  and  of  continuity.  This  is  in  contrast  with  the  more  common 
Eulerian  jellied- out  fluid  approach  where  the  particle  identity  is  lost.  The 
model  used  is  developed  in  some  detail.  Trajectory  plots  are  given  for  the 
faster  and  slower  space-charge  waves  by  themselves  and  for  velocity 


modulation  for  the  ratios  of  driving  frequency  to  plasma  frequency, 
0.  5,  1.  0,  and  2.  0. 
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FUNDAMENTAL  EQUATIONS,  MODEL 


In  one  formulation  of  the  Lagrangian  form  of  hydrodynamic  equa¬ 
tions,  the  independent  variables  are  the  time,  t,  and  the  initial  (t  =  0) 
coordinates  of  any  particle  of  the  fluid,  say  a,  b,  c.  The  dependent 
variables  are  the  particle  coordinates,  say  x,  y,  z;  their  values  in  terms 
of  a,  b,  c,  t  give  the  whole  history  of  every  particle  of  the  fluid.  The 
components  of  velocity  and  acceleration  of  a  given  particle  are  given  by: 


8t 

respectively 


8x  8  y  8 z 
'  8  t  ’  8  t 


and 


8^x 


8t^ 


8  2z 


8t“  8 8t 

Hence,  the  x  component  of  the  equation  of  motion  is 


m 


8  2x 
8  t^ 


=  F 


(1) 


where  F  is  given  at  the  position  of  the  particle,  m  is  the  particle  mass. 
The  equation  of  continuity  is  not  the  same  as  that  in  the  Eulerian  formulation. 

The  quantities  a,  b,  c  need  not  be  the  initial  coordinates  of  a  par¬ 
ticle,  but  may  be  any  three  quantities  which  serve  to  identify  a  particle  and 
which  vary  continuously  from  one  particle  to  another.  In  the  one-dimen¬ 
sional  model  to  be  used  here  it  is  desired  to  obtain  the  trajectories  of 
particles  beyond  the  plane  of  injection  say,  x  =  0.  Each  particle  will  be 
identified  by  its  time  of  departure  from  x  =  0,  the  a^^  particle  by  t  , 
etc.  Thus  for  the  a  particle  x  =  x(t,  t  )  The  result  is  to  be  a  plot  of 
particle  position  as  a  function  of  time  t,  such  as  shown  in  Figure  1. 

The  notion  of  a  particle  per  se  in  a  continuous  fluid  may  present 
some  difficulty  in  conception.  Of  course,  in  a  fluid  of  electrons  and  ions, 
one  might  be  tempted  to  ask  for  the  history  of  each  charged  particle; 
however,  this  request  is  absurd  as  there  may  easily  be  10^^  particles 
in  a  given  problem.  An  alternative  is  to  follow  selected  elements.  The 
selection  must  be  made  with  some  care  in  order  to  obtain  correct  results. 

In  the  continuous  (Eulerian)  fluid  model,  the  dependent  variables, 

V,  p,  i,  describing  the  electron  motion,  can  be  shown  to  vary  with  time 
as  cos  w^t  as  seen  by  an  observer  riding  with  the  stream  at  average 


-1- 


velocity  v^  -  The  implication  is  that  each  element  of  the  fluid  is  performing 
simple  harmonic  motion  at  the  radian  frequency  independent  of  the 
amplitude  of  excitation,  in  addition  to  its  drift  at  velocity  v^>  Thus,  a  model 
of  selected  elements  or  discrete  particles,  which  is  chosen  to  approximate 
the  continuous  fluid  model,  should  show  similar  motion. 

The  one -dimensional  model  to  be  used  here,  consists  of  electrons  in 
a  neutralising  background  of  positive  charges.  One  choice  is  to  let  the  fluid  of 
electrons  and  the  ions  be  represented  as  a  manageable  number  of  discrete 
particles,  that  is,  as  negatively  and  positively  charged  sheets.  The  sheets  may 
be  thought  of  as  locations  of  point  charges  of  very  low  density  and,  therefore, 
may  oscillate  through  one  another.  Starting  from  near  uniformity  (i.  e. ,  + 

and  -  sheets  interlaced  and  equally  spaced  from  each  other),  with  initial 
displacements  less  than  the  sheet  spacing  there  will  be  no  motion;  this  is 
unlike  a  plasma.  Starting  with  initial  velocities  just  large  enough  to  cause 
nearest-neighbor-crossovers,  the  motion  will  be  oscillatory,  to  be  sure, 
but  is  not  simple  harmonic  (not  sinusoidal,  but  anharmonic)  and  the  period 
is  dependent  on  the  amplitude  of  excitation.  In  order  to  make  the  electron 
sheets  perform  simple  harmonic  motion,  as  in  the  fluid  noodel,  each  electron 
sheet  must  be  given  sufficiently  large  excitation  to  cross  over  many  ion  sheets 
in  a  period-  For  trajectory  calculations,  accounting  for  the  crossovers  is  a 
bit  messy  because  the  force  jumps  at  each  crossover.  However,  note  that 
real  ions  are  much  heavier  than  real  electrons  (m^/m^  ^  1836)  so  that 
ion  motion  is  generally  very  small  compared  with  electron  motion,  so 
small  that  it  will  be  neglected  here.  Hence,  the  ion  motion  is  "solved" 
and  the  ions  no  longer  need  to  be  made  discrete  but  can  be  left  as  a  uniform 
fluid  or  background.  In  this  model,  the  electron  sheets,  even  for  the 
smallest  excitation  will  cross  over  ions,  and  execute  simple  harmonic 
motion.  The  remaining  restriction  in  the  present  calculation  is  that 
the  electron  sheets  are  driven  weakly  so  that  they  do  not  cross  each 
other,  which  is  equivalent  to  the  requirement  of  a  single-valued 
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electron  velocity  in  the  fluid  model.  (Large  drive  will  cause  electron 
crossovers,  making  a  direct  solution  difficult,  but  easily  handled  on  a 
high-speed  computer.  ) 

The  model  to  be  analyzed  is  sketched  in  Figure  2.  A  portion  of 
a  typical  charge  density  and  electric  field  pattern  is  shown  in  Figure  3, 
assuming  electron  sheets  of  finite  thickness. 

The  detailed  motion  is  found  readily  from  the  isolated  electron- sheet, 
ion-block  pair  shown  in  Figure  4.  The  total  charge  is  assumed  to  be  zero, 
so  that 

■Pe^^e  =  Pi^^i 


The  ion  block  has  been  assumed  to  be  incompressible;  d^  is  constant.  The 
electron  sheet  also  is  considered  incompressible  so  that  it  moves  under 
the  influence  of  the  force  spatially  averaged  over  the  sheet.  The  center  of 
the  sheet  sees  average  field  IT  as  shown,  given  by 


E  = 


(X  .  x^) 

o 


d.  +  d^ 
__i _ e 

2 


(3) 


Because  the  field,  hence  force,  is  proportional  to  displacement,  it  is 
derivable  from  a  potential  which  is  parabolic,  V««*(x  -  x^)  .  In  the  equation 
of  motion,  the  charge-to-mass  ratio  of  the  sheet  appears;  this  is  simply 
e/m  of  a  single  electron.  Thus,  this  equation  is 


(X  -  x^) 


(X 


*o> 


(4) 


This  is  the  equation  of  simple  harmonic  motion,  as  desired.  The  plasma 
frequency  used  here  is  defined  by 


(5) 


using  the  electron  charge  and  mass  but  the  ion  background  density;  the 
sheet  and  block  thicknesses  do  not  appear  explicitly  and  may  be  made 
arbitrarily  large  or  small,  e.  g.  ,  d^  — >  0 

Lastly,  many  pairs  nnay  be  placed  one  after  the  other,  making  the 
background  continuous  and  uniform  for  equal  p.'s  The  x^'s  will  be  taken 
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ELECTRON  SHEETS 
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X, CENTER  OF  SHEET 

Figure  4 
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arbitrarily  as  the  center  of  each  block,  an  obvious  equilibrium  point,  so 
that  the  above  equation  of  motion  still  holds.  However,  with  the  background 
now  continuous,  any  sheet  may  leave  its  parent  block;  the  only  restriction 
imposed  in  this  study  is  that  sheets  may  not  cross  each  other.  Thus,  the 
notion  of  blocks  may  be  set  aside,  and  one  need  keep  track  of  only  the  sheet 
and  the  equilibrixim  positions. 

The  block-sheet  model  may  be  extended  to  other  symmetries. 
Starting  with  a  uniform  quiescent  sea  of  electrons  and  (immobile)  ions  and 
nnaking  a  local  displacement  of  electrons  along  one  coordinate  r,  one  can 
readily  obtain  the  equation  of  motion.  Let  the  equilibrium  radius  be  r^  and 
the  displaced  electron  shell  be  at  r  =  r^  (1  +  &)■  Then  the  electron  motion  is 
governed  by 


I?"  n+ 1  L  (1  +  J 


n  =  0,  1,  2  for  planar,  cylindrical  and  spherical  displacements,  respectively. 
The  motion  will  be  oscillatory  but  not,  in  general,  simple  harmonic.  For  all 
planar  displacements  (n  =  0)  but  only  for  small  cylindrical  and  spherical 
displacements  (n  =  1,2,  6  « 1)  this  equation  reduces  to 


(7) 


Thus,  simple  harmonic  motion  is  obtained  only  under  these  restrictions. 
See  Dawson  (1959)  and  Jackson  (1960). 


PARTICLE  DISPLACEMENT 

Particle  orbits  are  obtained  by  integrating  the  equation  of  motion 
of  each  particle.  The  general  answer,  for  the  a^^  particle  is 

x=  x{t,  t^)  (1) 

=  +  F(t^)  cos  «p(t  -  t^)  +  G(t^)  sin  «p(t  -  t^)  (2) 
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The  F's  and  G's  are  the  two  independent  constants  allowed  by  the 

second  order  equation  of  motion;  they  may  be  arbitrary  functions  of  t  • 

th  th  ^ 

A  similar  equation  will  hold  for  the  b  ,  c  .  .  .  particles,  each  with  its 
own  x^,  F  and  G.  The  displacement  of  each  particle  is  independent  of 
the  others  so  long  as  there  are  no  crossings,  i.  e.  ,  the  order  of  the  sheets 
is  maintained.  If  two  sheets  do  cross,  they  must  exchange  x  's,  making  a 
jump  in  the  force;  velocity  and  displacement  are  continuous  through  the 
crossover. 


DRIFT  VELOCITY 


Each  particle  executes  simple  harmonic  motion  about  its  equilibrium 
position,  x^.  It  is  desirable  to  obtain  the  motion  of  particles  in  a  drifting 
stream  from  these  results.  Hence,  the  equilibrium  positions  will  be  made 
to  drift  at  velocity  v^,  that  is. 


X  =  V  (t 

o  o 


‘oa) 


(1) 


Or,  the  original  frame  of  observation  may  be  considered  to  be  drifting  at 

velocity  v^  so  that  a  simple  (Galilean)  transformation  may  be  made  to  the 

laboratory  frame  to  obtain  the  drift-plus-oscillatory  motion,  t  is  the 

th 

departure  time  of  the  observer's  frame  for  the  a  electron  and  is 
different  from  t^  if  there  is  an  initial  displacement.  On  occasion,  it  may 
be  inconvenient  to  carry  t^,  a  fictitious  departure  time  and  hence,  it  is 
desirable  to  have  a  way  to  eliminate  t^.  Using  x  =  0  at  ^ 
displacement  equation  is 


»  =  ''o'*.  -  'o.'  * 
From  this  result,  one  forms 


x(t,  t^)=  v^(t  -  t^)  +  F(t^)  |cos  «p(t  -  t^)  -  ij  +  G(t^)|8in  w^Ct  -  t^)j 
The  pertinent  times  are  shown  in  Figure  5. 


(2) 

(3) 
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EQUILIBRIUM 
POSITION  ■ 


TRAJECTORY- 


F(ta)co8«p<t-ta) 
•»0(to  )8in  ft»p(t-ta) 


F(tfl) 


•Vo(t“t0(|) 


One  might  aik  whether  t^  and  vary  with  t^;  that  is,  may  the 
ion  blocks  depart  at  arbitrary  times  with  arbitrary  velocity?  In  general, 
the  answer  is  yes,  but  this  would  mean  letting  the  background  vary.  In  the 
restricted  uses  here,  the  answer  is  no.  The  individuality  of  the  blocks  is 
to  be  suppressed  and  the  background  charge  density  is  to  remain  a  constant. 
Thus  all  blocks  are  alike,  departing  uniformly  in  time  and  with  the  same  v^ 
However,  as  F(t^)  and  G(t^)  are  not  restricted,  one  is  still  quite  free  to 
put  in  arbitrary  excitation  of  the  electrons  including  steps,  ramps  etc. ,  as 
well  as  periodic  functions. 

INITIAL  VELOCITY 

The  constants  F  and  G  may  be  determined  in  several  ways.  Here 
the  choice  is  to  specify  two  independent  qualities  at  the  input  plane,  the  ac 
velocity  and  ac  current  density.  The  velocity  of  the  a^^  particle  is 

v(t,  =  =  v^  -  F(t^)«pSin  Wp(t- t^)  + G(t^)WpCOS  «p(t- t^)  (1) 

At  the  input,  x  =  0,  f  velocity  is  given  by  a  constant  value  plus  a 

time>dependent  part, 

(The  constant  value,  v^,  is  the  same  as  that  of  the  equilibrium  positions 
but  is  not  necessarily  the  average  value  in  the  presence  of  modulation.  ) 
Hence, 

Vi(ta)  =  G(ta)«p  (3) 

determining  one  of  the  constants.  Note  that  inay  be  quite  arbitrary, 

periodic  or  otherwise. 
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INITIAL  CURRENT  DENSTY 


The  current  density  or  particle  flux  depends  on  the  rate  of  passage  of 
electrons  through  an  area;  here  the  current  density  is  simply  proportional  to 
the  density  of  trajectories  at  a  fixed  x.  Similarly,  the  charge  density  is 
proportional  to  the  density  of  trajectories  at  a  fixed  t.  Hence,  if  the 
problems  "start"  at  a  fixed  plane,  as  occurs  here,  the  obvious  second 
initial  condition  to  be  specified  is  the  current;  however,  for  a  problem 
starting  at  a  given  time,  specifying  charge  density  would  be  more  fitting. 

The  calculation  current  density  can  be  naade  with  the  aid  of  Figure  6, 
showing  two  adjacent  trajectories.  The  di  rect  or  steady  current  density, 
i^,  is  inversely  proportional  to  t^  -  t^,  or 

,  _  U) 


p^  is  the  charge  /  area  of  the  electron  sheet,  i*  negative 

number.  The  total  current  density  depends  on  the  actual  time  between  sheets, 
which  is,  at  X  =  0,  using  Eq.  (2)  under  Drift  Velocity, 


As  F(t^)  is  a  continuous  function  in  the  limit,  it  maybe  expanded  as 

dr(t  ) 

F(t^)=  +  -  t^)+  -  -  -  (3) 


Putting  this  into  (2)  and  solving  for  t^  >  t^,  one  obtains 

^•••■'■•‘  ■-'[■•4;  ??]■' 


Hence,  the  current  density  at  x  =  0  is 


(4) 
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ta  too  ^ob 

Figure  6 
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(5) 


(6) 


At  the  input,  x  =  0,  ^  current  density  is  given  by  a  constant 

value  plus  a  time  -  dependent  part, 


‘a>  =  'o  W 


(7) 


(The  constant  value,  i^,  is  the  current  density  for  no  modulation  but  is 
not  necessarily  the  average  value  in  the  presence  of  modulation.  )  Thus, 
one  has. 


ii(v=io 


dF(t^) 

dt 


(8) 


Thus  (4)  becomes  simply. 


(9) 


Integrating  (8),  F(t  )  is  given  by, 
& 


(10) 


The  integral  is  made  definite  by  saying  that  for  t^  <  t^  there  is  no  current 
modulation,  and  hence,  no  displacement  from  uniformity;  the  current 
modulation,  periodic  or  otherwise  starts  at  t^  =  t^.  Thus  F(t^)  is  an 
accumulated  shift  from  no  modulation.  Physically,  i.dt  is  the  ac  charge 
density  added  or  taken  away  from  the  unmodulated  flow  and  the  integral 
is  the  accumulated  charge  added  or  taken  away.  In  the  usual  case  of  a 
periodic  input  current,  i|(t  1,  of  course,  maybe  given  by  a  Fourier  series, 

X  A 

with  a  typical  term  i  cos  nwt  ;  then  F'(t  )  must  have  the  same  term,  and 

n  a  a 

F(t  )  must  have  the  term 

ci 

on. 

I-  —  sin  n«t^ 


and  the  transient  associated  with  the  start  of  modulation  at  t^  generally 
can  be  neglected. 

Starting  from  fundamentals,  not  using  the  preceding  results,  one 
may  show  that  the  initial  ac  charge  density  is  given  by 

••o  r'*'*'.*  1 

<V  '»> "  7; 

The  current  at  any  plane,  i(t,  t^),  may  be  found  in  terms  of  F(t^), 
G(t^)  and  (t-t^).  The  interested  reader  is  advised  to  start  from  a  statement 
of  conservation  of  charge, 

Kta-  =  i(t.  1^)  dt. 


LIMITS  ON  EXCITATION 


In  the  cases  of  interest  here,  the  excitation  is  to  be  maintained 
sufficiently  small  so  that  sheets  do  not  cross  over,  that  is,  the  ordering 
of  sheets  is  to  remain  fixed.  At  a  fixed  x,  going  to  larger  t,  one  must  find 
sheets  that  departed  later,  that  is 


(1) 


This  ratio  may  be  calculated  from 


8t 


X 


_  ax/st 

TITir 


Or,  at  a  fixed  t,  going  to  larger  x,  one  must  find  sheets  that  left  earlier, 
that  is 


(2) 
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Or,  ,in  a  problem,  which  starts  at  a  fixed  time,  t  =  0  with  electrons 
initially  ordered  as  x^,  x^  etc.  t  ordering  is  maintained  at  a  later  time 
if 

>  0  (3) 

t 

(In  this  problem  t  -  t  is  replaced  by  t  -  0  =  t  and  F(t  ),  G(t  )  are  replaced 
by  F(x  ),  G(x  )■  At  crossings,  the  densities  become  infinite. 

As  brought  out  in  the  examples  to  follow,  the  excitations  can  be  quite 
large  before  overtaking  occurs,  considerably  larger  than  allowed  in  the 
small-signal  Eulerian  fluid  analysis  where  higher  order  terms  are  dropped, 
and  no  harmonics  are  allowed.  Here,  no  such  restrictions  occur.  The  only 
problem  is  with  overtaking  and  that  could  be  obviated  by  renumbering. 

Example  (!)■  Assuming  a  velocity  modulation  of 

'^l<‘a>  =  '^1  “‘a 

crossing  just  occurs  for  u>w  with 


between  the  electrons  that  left  in  the  vicinity  of  -  ir/ Z  (or  3ir/2,  etc.) 

at  u  (t-t  )  =  ir/ 2,  one  quarter  plasma  cycle  after  modulation.  For  u<u  , 

P  ®  P 

v^/v^  >  1  causes  crossings  at  the  input  between  electrons  near  ut^  =  -  ir. 

Thus,  one  may  plot  the  nuiximum  allowable  velocity  modulation  as  in 

Figure  7a. 

Example  (2)-  Assuming  a  current  modulation  of 

u 

crossing  just  occurs  for  a*  >— £  at 

2 


=  0.  5 


between  the  electrons  leaving  near  (wt^  *4)^  ir  at  the  input.  For 
ci)<***p  crossing  just  occurs  for 
2 
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Figure  7a 


for  crossings  between  (cat  -  ^  +  ir/2  electrons  at  c»  (t  -  t  )  =  ir/2  one 

p  A 

quarter  plasma  cycle  after  modulation.  These  results  are  shown  in 
Figure  7b. 

Example  (3).  For  excitation  of  the  slower  space -charge  wave  only,  the  limits 
of  excitation  in  velocity  and  current  are 


P  ft  ~  '  D 


as  shown  in  Figure  7c. 

Example  (4).  For  excitation  of  the  faster-wave  only  for  ca<  ca^ 


and  for  u>>ui 

P 


^<l/(-l+  )  /  ( -1+ 

o  P  ^O  P 


as  shown  in  Figure  7d. 


RELATION  TO  FLUID  MODEL. 


The  particle  results  may  be  put  in  the  form  that  goes  with  the 

continuous  fluid  model  used  in  the  original  space -charge  wave  work  of 

Hahn  and  Ramo  (1939).  In  the  present  form,  one  has  the  particle  position 

x(t,  t^):  in  the  fluid  model,  the  dependent  variables  are  given,  for 

example,  as  v(x,  t)-  The  transformation  from  one  system  to  the  other  is 

obtained  by  going  to  a  frame  moving  at  velocity  v^,  as  was  mentioned 

earlier,  and  observing  the  streaming. 

An  observer  leaves  x  =  0  at  t  =  t  with  velocity  v  and  has 

oa  '  o 

position 

*  =  -  ‘oa> 
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(1) 
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Figure  7c 
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Figure  7d 


(2) 


Or,  he  measuree  time  by 


t  =  ~ 

oa  V 

o 

X  +  F  (t  ) 

=  t  + - L. 

a  V 


(3) 


The  independent  variable  t^  may  now  be  replaced  in  terms  of  t  and  x 
in  the  particle  results, 


r  X  +  F  (t  )  -| 


(4) 


or 


sin  u  (t  -  t  ) 
P  a 


(5) 


However,  the  trigonometric  terms  multiply  ac  (first-order)  terms,  so 
that  in  the  limit  of  very  small  current  modulation,  the  ac  term  F(t  ) 
in  the  arguments  may  be  dropped.  Hence,  the  ac  displacement,  as  seen  by 
the  observer,  is 


Xj  =  Xj  (x,  t) 


(6) 


W  x  W  X 

F  (t  -  ~  )  cos  -J2-  +  G  (t  -  —  )  sin  -E- 


(7) 


The  wavelike  nature,  which  was  completely  suppressed  in  the  particle 
picture,  is  once  more  evident.  However,  here  one  treats  x  and  t  as 
independent  variables,  no  longer  connected  together,  as  would  be  so 
for  the  trajectory  of  a  given  particle.  The  ac  velocity,  as  seen  by  the 
observer,  is  obtained  from 


Vj  =  Vj  (x,  t) 

*^1  . 

=  •ar 

8x 


— L  +  V 
at  o  ax 


(8) 

(9) 

(10) 
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F  and  G  are  functione  of  a  single  variable  so  that  one  may  write 


\  ^  -  pt  pi 

TT-  TT-  -T  ^ 

o 

Thus 

Vj  (X,  t)  =  «  G  (t  -  i  )  cos  „  p  (t  .  i  )  ,in 

O  0*^00 

At  X  =  0,  let 

Vj  (0,  t)  =  cos  w  t 

=  o  l«  -  ) 

or 

G  -  ^  )  =  CO*  ^  ) 

o  o 


(11) 

(12) 

(13) 

(14) 

(15) 


The  current  density  i,  (x,  t)  must  obey 

^  dE, 

i^(x,  t)  =  - 


(16) 


was  obtained  earlier  as 


o 


*1 


(17) 


(The  ion  charge  density,  p^,  is  the  negative  of  the  electron  charge 
density,  p^,  in  the  fluid  model. )  Hence, 
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(18) 


ij  (x. 


t)  = 


P 


o  8  t 


At 


= "o  r  (‘  -  f  )“• 

fa>  X 

<i>  X  n 

^  +  G'  (t  -  -^  ) 

(19) 

L  o 

o  o 

O  J 

0,  let 

il(0.  t) 

=  ij  cos  <(i>t  -  4) 

(20) 

=  p„  (t  -  f  > 

(21) 

o 


or,  integrating, 

F  (t  -  ~  ) 


•in  (ut 


(22) 


These  results  may  be  summarized,  as 


V 


V 


(i>  X  1, 

.-E_  1 

V  T 


W  X 

-E_ 

V,  V  U  X 

+  —  COS  (fa)t  -  — )  sin  — ^ 

V 

Q 

< 

> 

c 

3 

c 

> 

V 

JT 

(23) 

h> 

^  (A)  X 

-E 

b> 

sin  (wt  -  —  -  ♦  )  •in  — 

o  '^o 

(24) 

ii(x>t)  i  w  X  V,  w  X 

- -  =  ^  cos(«t-—  -♦)cos^.^  ^sinM-^  )  sin-^  (25) 

i  i  o  o  o  p  o  o 


p.  may  be  obtained  from 

^  fb  C 


Pi  (*• 

‘o“yir 

8x 

-  ■  Po“53r 

(26) 

(27) 

The  result  is 

Po 

r 

W  X 

u  x“] 

Pi  (X.  t)  =  :^ 
o 

I  (F'  -  G»p )  cos 

^o 

+  (F  »  +  G')  sin  — 

P  Vo  J 

(28) 
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Hence 


(w  t  -  —  ) 
o 


cos 


W  X 


o 


sin 


sin 


W  X 


o 


(29) 


These  results  are  the  same  as  obtained  for  the  fluid  model. 

One  final  word  of  caution;  x^  (x,  t)  as  given  here  is  the  displacement 
from  equilibrium  of  electrons  as  seen  by  an  observer  at  (x,  t)  and  is  not  the 
displacement  of  a  particular  electron. 


TRAJECTORY  PLOTS 


The  properties  of  various  kinds  of  modulation  and  of  each  wave 
may  be  made  clearer  using  plots  of  electron  trajectories.  The  method  of 
obtaining  these  plots  will  be  given  in  some  detail  for  harmonic  excitation. 
The  starting  point  is  with  the  displacement, 

x(t.  t^)  =  %5(t.  t^)+  F<t^)  l^cos  «p(t-  t^)j+  G(t^)  l^sin  «p(t-  t^)j 

The  initial  conditions  are  stated  as 

v^  (t^)  =  v^  cos  wt^  ( 1 ) 

ij  (ta)  =  i^  co«  +  )  (2) 


v^  is  real  and  positive  and  i^  and  i^  are  real  and  negative.  The 
displacement  is 

i,  V 

x(t,  t)=v(t-t  )  +  -r-—  sin  ( wt  -  f )  cos  w  (t  -  t  ) 
^  a'  o  oa  1  u  a  ^  p  a 

o 

V,  V 

+  7“  IT  "‘a  "p  ■  ‘a> 

o  p 


(3) 
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Or  the  x^(x,  t)  expreeaion  may  be  used  with  x  and  t  properly  related  in  the 
functions  so  as  to  follow  a  particular  electron,  i.  e. , 


X  =  V  (t  -  t  .) 

o  oa' 


For  plotting,  one  needs  to  have  a  schedule  of  departure  times.  These 
times  occur  uniformly  in  t^,  but  not  in  t^.  The  electron  departure  times 
are 


=  'o.  -  4-  '■'‘.I 

O 

1 

=  t  -  .ri  i.  sin  (wt  -  4) 
oa  1  o>  '  a  ^ ' 

o 


(4) 

(5) 


This  is  a  transcendental  equation  in  t^  which  can  be  solved  numerically: 
for  small  current  modulation,  it  can  be  solved  approximately  as. 


t  »t 
a  oa 


ji  i  .in  (o.t^  -♦) 
o 


(6) 


The  t  's  may  be  scheduled  as 
oa  ' 

t  =  t  =  nAt 
oa  on 


n  —  1,  2,  3  .  •  • 


(7) 


where  At  is  the  uniform  spacing,  at  a  given  x,  of  the  ion  blocks  or 
equilibrium  positions.  For  N  electrons  injected  in  each  cycle,  ut=2ir, 


N  At  =  2  v 


(8) 


so 


t  =  2«-^ 
on  « 


(9) 


Hence, 


X  (t,  X  )  — e  X  (t,  n)  or  x  (t) 
a  n 


(10) 


For  convenience,  the  passage  of  time,  ca  (t>  t^),  may  be  measured 


in  terms  of  an  equivalent  distance.  Let  x^^  be  defined  by 


X  =  V  (t*t  ) 

on  o  '  a' 


n^**  electron 


(11) 
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(Thia  is  the  equilibrium  poeition  only  in  the  case  of  no  current 
modulation. )  Then, 


“p*on 

«o(t-t  )  = - 

pa  V 


(12) 


Hence  x(t,  t^)  becomes 


WpX 

v_ 


(x^^.  t^)  =  Wp(t  -  t^^)  +  ^  ^  sin  («„t^  -  cos  ( 


ti>  X 

p  on 


+  -  cos  «t  sin  {-£ - 

V  a  V 

o  o 


) 


(13) 


In  this  form,  trajectory  plotting,  x  vs  t,  is  straightforward: 

By  way  of  illustration,  one  may  find  the  locus  of  turn-around 
points  (where  the  ac  velocity  changes  sign).  For  velocity  modulation 
only,  t^  =  t^,  and  particles  obey 


V.  V 

x(t,  t^)  =  v^(t- t^)  + - 2-  cos  wt^  sin  w  (t-t^) 

o  p 

v(t,  t^)  =  v^  +  Vj  cos  wt^  cos  w  (t  -  t^) 


(14) 


(15) 


The  maximum  excursion  from  equilibrium  occurs  for  u>p(t  -  t^)  s  w/2 
who  re  the  ac  velocity  vanishes,  v  is  v^;  but  the  electron  at  turn  around 
is  not  at 


w  X 

P 


of  v/2,  rather,  at 


w  X 

P 


turn-around  =  —  + 
2 


I 


cos  wt 


(16) 


'  o  '  o 

This  locus  and  some  representative  trajectories  are  shown  in  Figure  8. 


SLOWER  SPACE-CHARGE  WAVE  TRAJECTORIES 


The  slower  space-charge  wave  by  itself  may  be  excited  by  proper 
choice  of  initial  modulation.  As  known  from  space- charge  wave  analysis, 
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this  requires  that  the  initial  velocity  and  current  modulation  be  ir  out  of 
phase,  ^  =  IT,  and  in  the  ratio 


Hence, 


(1) 


u 


•,)  =  "p  <‘  -  ‘o.> 

O 


=  u  (t  -•  t  ) 
p  '  oa' 


'^l 

V 

o 

sin 

«ta  - 

P  a 

sin 

‘a 

b>  X  "] 

P  on 

V 

o 

"o  J 

(2) 


Note  that  all  all  electrons  have  the  same  amplitude  of  oscillation. 

Figures  9,  10,  and  11  are  plots  of  trajectories  of  the  slower  wave 

for  w/Up  of  2,  1,  0.  5.  The  regions  of  larger  charge  density  are  where  the 

electrons  are  going  slower,  giving  rise  to  negative  ac  kinetic  energy  density. 

The  phase  velocities  for  the  figures  are  2/3  v  ,  1/2  v  and  1/3  v  , 

o  o  o 

obviously  corresponding  to  the  velocities  of  propagation  of  the  denser  regions. 
The  bunches  are  formed  then  reformed  with  electrons  sliding  through  the 
bunches.  There  is  no  collection  of  particles  which  stays  more  or  less 
together;  that  is,  there  is  no  trapping  of  particles. 


FASTER  SPACE- CHARGE  WAVE  TRAJECTORIES 


Excitation  of  the  faster  space-charge  wave  by  itself  requires  that  the 
current  and  velocity  modulation  be  in  phase,  ^  =  0  and  in  the  ratio 


V 

o 


!!!e 

w 


h 

^o 


(1) 


Hence, 
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^  «(t,  t  )=»  + 

o 


'  «p<‘  -  ‘o»)  ♦ 


sin  +  «p  (t  ’  J 

i  •“[■■•  •=?] 


Again,  the  amplitudes  of  oscillation  are  all  the  same.  Figures  12,  13,  14 
show  the  trajectories  for  =  2,  1,  0.  5.  Here  the  regions  of  larger  charge 

density  occur  where  electrons  are  going  faster,  giving  rise  to  positive  kinetic 
energy  density.  The  phase  velocities  are  2v^,  oo,  and  -v  ,  not  so  obviously 
tied  to  the  velocities  of  the  denser  regions  as  is  true  with  the  slower  wave. 

For  Up,  only  i^^  and  v^£  are  excited  and  is  zero;  hence  at  fixed  t, 
the  spacings  of  the  trajectories  remain  fixed.  In  addition,  as  p^^  =  0, 

•  E/tx  =  0,  so  Ej^  =  E^(t)  independent  of  x;  thus  the  electric  Held  over 
all  space  is  merely  oscillating  with  time  at  u  =  u  . 


VELOCITT  MODUUITION,  WAVE  INTERFERENCE 


For  velocity  modulation  only,  both  waves  are  excited.  The  displace¬ 


ment  is. 


(t,  t  )  =  «  (t  -  t  )  +  —  cos  b>  t  sin  fa>  (t  -  t  ) 
a  p  oa  V  a  p  '  a 

r  Q 


Plots  of  these  trajectories  are  given  in  Figures  15,  16,  17  for  w/Up  of  2, 

1.  0,  0.  5.  Here  the  wave  interference  between  the  faster  and  slower  waves 
is  evident.  The  slower-wave  tends  to  be  predominant,  due  to  absence  of 
p^  at  w  =  Up  and,  in  the  other  two  plots,  because  =3  p^_  . 

Vj^  and  v^  are  equally  excited  as  are  and  i^  at  all  frequencies. 

Although  the  current  density  is  always  a  maximum  at  from 

the  plane  of  excitation,  independent  of  u/Up,  the  charge  density  follows  the 
sine  law  only  for  large  u/Up.  In  the  small  WWp  region,  the  charge  density 
follows  a  cosine  law  and  gives  the  appearance  of  strung  out  bunching,  as 
seen  in  the  figures. 
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CUULENT  MODULATION 


Zb  the  c&ee  of  current  modulation  only, -the  displacement  is 


(iX  ^ 

*(t,  t^)  =  w  (t  -  t^)  +  -p  -jf  -  ♦)  cos  «  (t  -  t^) 


(1) 


This  may  be  compared  with  the  case  for  velocity  modulation.  Let  the 
preceding  t^  time  axis  be  shifted  so  that  (wt^  •  4)  becomes  (cat^'-t- v/2). 
Let  the  (t-t^)  axis  be  shifted  so  that  (t  -  t^)1  becomes 


The  new  displacement  is 


w  X 


i,  «_ 


V 

o 


( t' .  t;  )  =  «  (f  -  f^)  -  */2  +  ^  cos  «;  sin  «  (f  - 1;) 


(2) 


p  oa'  “  P 

This  is  the  same  form  as  for  velocity  modulation,  but  shifted  in  time  and 
displaced  v/ 
unnecessary. 


displaced  v/2  in  normalised  distance  {k^/4).  Hence,  new  plots  are 


INTERACTION  WITH  A  TRAVELING  WAVE 


The  stream  may  react  to  its  own  fields,  space-charge  fields,  and 
may  also  be  driven  by  externally  applied  fields.  Such  fields,  of  course,  are 
necessary  to  provide  the  initial  modulation,  v^  and  i^,  used  earlier; 
these  fields  presumably  are  imposed  over  a  short  length  of  stream,  a  diode 
or  gap  region.  Particle  motion  in  a  gap  has  been  given  in  great  detail  by 
many  authors  over  the  past  three  decades,  as  found  in  the  classical 
monograph  by  Llewellyn  (1941)  and  as  outlined  by  one  of  the  earliest 
workers  with  transit  time  effects,  Benham,  in  Benham  and  Harris  (1957). 

Interaction  with  a  traveling  wave  provides  an  interesting  example 
of  a  driven  stream,  as  well  as  insight  into  the  initial  bunching  process  in  a 
traveling- wave  tube.  The  problem  solved  here  is  very  simple  and  somewhat 
non-physical.  A  slow  wave  of  phase  velocity  v^  is  caused  to  travel  with  a 
stream  of  average  velocity  v^.  The  wave  produces  bunching  and  may  alter 
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the  stream  kinetic  energy:  this  change  in  kinetic  energy  will  appear  in  field 
energy,  some  of  which  may  cause  the  driving  wave  to  increase  (as  in  a 
traveling- wave  tube)  or  decrease  (as  in  an  accelerator),  all  as  a  self-con¬ 
sistent  process.  However,  here,  only  the  initiation  of  bunching  will  be 
obtained  because  the  driving  field  will  be  taken  to  be  constant. 

The  electric  field  of  the  driving  wave  is  given  in  terms  of  laboratory 
coordinates  as, 

EJx.  t)  =  sin  («  t  -  ^  )  (1) 

P 


The  equation  of  motion  for  a  particle  is 


m  0^x 


at* 


e  E  .  +  e  E 

space  charge  wave 


(2) 


where  the  fields  are  given  at  the  particle.  The  particle  position  x  appears 

in  E  ,  whereas  solution  would  be  easier  if  the  equilibrium  position  x 
wave  ^  o 

appeared.  Hence,  one  may  expand  E^  as 


8E  (x  ,t) 

=w  (X. (x„,  0  t  ,  (*..„)+--- 

=  Ew  -  7“^  )  -f-  )  E  cos  ( wt  -  ^  )  + 

P  P  P 


The  phase  is 


p  w 


(3) 

(4) 

(5) 


where  is  the  wavelength  of  the  driving  wave.  In  the  limit  of  small 

excitation,  the  deviation  x-x  will  be  much  smaller  than  X  ,  so  that 

o  w 

the  second  term  may  be  ignored.  The  problem  solved  is  for  no  initial 
current  modulation  and  no  initial  velocity  modulation.  Thus,  t^  =  t^. 
Hence 

*o  =  '^o  ■  ‘a) 


-41- 


(7) 


Thu*,  the  phase  of  E  becomes 

p  p 

=  «(i  +  ^  («) 

p  p 

V  (lit 

=  M't+  -SL-*  (9) 

'^P 

w'  is  the  frequency  of  the  wave  as  seen  by  the  stream,  the  usual  Doppler 
shifted  value.  The  response  of  the  stream  to  the  wave  is  expected  to  be 
greatest  if  ta*  is  at  the  resonant  frequencies  of  the  stream,  t  or 

(1  -  ^  )  =  t  «  (10) 

^  ^  ^  (11) 

7“  ''  V  V 

poo 


That  is,  greatest  response  is  expected  if  the  wave  velocity  is  that  of  the 
faster  or  the  slower  space-charge  wave.  Physically  this  requires  that  the 
wave  should  move  one  wavelength  past  the  stream  in  one  cycle  of  the 
plasma  frequency:  this  slip  is  just  that  used  to  obtain  maximum  growth 
rate  in  traveling- wave  tubes  with  mode  rate -to- large  space-charge  den-  . 
sities.  Choosing  this  synchronism,  one  obtains  the  equation  of  motion. 


"p  =-^=w  -  "p  ‘a)  +  "‘aj 


The  solution  is 


*(t,  t^)  = 


COS  wt 
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The  (A>  term  shows  the  growth  of  displacement  from  equilibrium  due 

to  traveling- wave  interaction,  increasing  linearly  with  time  of  flight;  however, 

the  motion  is  still  oscillatory,  with  electrons  sliding  through  the  bunches  and 

not  being  trapped.  If  the  synchronism  is  with  the  slower  space- charge  wave 

(lower  sign),  then,  bunches  will  be  formed  in  the  retarding  phase  of 

this  implies  decrease  of  kinetic  energy  of  the  stream  and  increase  of  wave 

energy,  that  is,  an  increase  in  E^-  Had  an  energy  balance  been  required,  E^ 

would  grow.  Then,  in  turn  would  grow  even  more  rapidly,  as  would  E^ 

resulting  in  the  usual  exponential  growth.  Similarly,  for  synchronism  with 

the  faster  space-charge  wave,  bunches  are  formed  in  the  accelerating  phase 

and  E  must  decrease,  as  in  an  accelerator, 
w 

Solutions  of  the  self-consistent  problem  would  be  more  useful.  Also,  in 
either  solution,  as  the  deflection  grows,  particles  must  cross  at  some  point  so 
that  complete  solution  to  maximum  growth  of  the  wave  or  stream  energy  virtu¬ 
ally  requires  use  of  a  high  speed  computer;  typical  nuichine  solutions  and 
trajectories  are  given  by  Tien  (1955). 

HISTORICAL  NOTE;  ACKNOWLEDGMENT 

Distance-time  plots  are  quite  commonly  made  for  one -dimensional 
problems.  Perhaps  the  most  well-known  electron  stream  plot  is  that  for  the 
klystron  as  drawn  by  a  patent  lawyer,  Applegate  (1942),  for  the  early  Varian- 
Hansen  klystron  patent  applications.  In  his  drawing,  space-charge  forces  were 
neglected.  There  is  the  temptation  to  use  the  Eulerian  space-charge  wave 
solutions  as  a  guide  to  altering  the  Applegate  diagram  so  as  to  show  the  effect 
of  space-charge  forces.  These  efforts  generally  fail  with  the  common  and 
easily-detected  error  of  putting  the  turnaround  points  at  ^pf^'  particle 

trajectories  shown  here  were  initially  solved  for  and  sketched  by  J.  L.  Palmer 
for  a  graduate  seminar  given  by  the  writer  at  the  University  of  California, 
Berkeley,  in  1959.  A  proper  name  for  them  would  be  Palmer  diagrams.  The 
derivation  given  here  is  somewhat  different,  following  the  classical  Lagran- 
gian  form  as  used  by  Benham  (1928)  and  Llewellyn  (1941)  for  diode  problems. 
However,  many  steps  here  are  parallel  to  those  of  Palmer  to  whom  the  writer 
is  most  grateful.  The  draivings  shown  are  due  to  the  computing  of  R.  Lanstein, 
the  careful  plotting  by  Betty  Jean  Hansen  and  the  skilled  drawing  of  George 
Faraco. 
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